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We discuss the theme on translating different descriptions of computative algorithms 
into high-level programming languages, enumerate some advantages of analytical 
descriptions and demonstrate that logical functions may be used effectively to create 
analytical formulae available for describing a set of combinatorial and number-theo- 
retic computative algorithms. In particular, we adduce analytical formulae to 
generate /-th prime numbers p,, permutations of order m, k-th numbers of 
Smarandache sequences of Ist and 2nd kinds and classical Magic squares of an 
order n. 
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1 Introduction 


As well-known"? verbal and diagram (graph-diagram) techniques available for 
describing computative algorithms are the most wide-spread at present. For instance, 
Euclidean algorithm, allowing to find the greatest common divisor (GCD) of the 
positive integers a and b (a > b) has the following verbal description ° 

1. Assign m=a,n=b; 

2. Find r=m mod n; 

3. Ifr > 0, then pass to 4. Otherwise, pass to 5; 

4. Assign m =n, n =r and pass to 2; 

5. Answer: GCD(a, b) =n. 

Since all computative algorithms are realised, as the rule, on computer at present, 
the main fault of the verbal description of computative algonthms is the necessity of 
translating this description into one of special computer-oriented languages. 

The diagram form of the description of computative algorithms allows to 
simplify slightly the process of such translation. In particular, the diagram form of 
Euclidean algonthm is shown in figure 1, where squares with digits 1, 2, 4, 5 and 
the rhomb with the condition r > 0 mean respectively to points 1, 2, 4, 5 and 3 of 
the verbal description. 


Figure 1. Diagram form of the description of Euclidean algorithm. 
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The logical technique‘ available for describing computative algorithms is less 
known than verbal and diagram descriptions, but just it gives easier way to obtain a 
program code. In particular, we may present the logical description of Euclidean 
algorithm in the form 


Ay7A203T AgBT2t As, (1) 


where A; are elements of the vector operator A = {(m = a, n = b), (r = m mod n), 
(@), (m=n,n =P), (n)}; a; are elements of the vector conditional jump @ = {(©), 
(B), (r>0), (BD), (D)}; (B) is the blank in A and a; 1, and 1’ are arrows 
indicating respectively points of departures and destinations; B is the unconditional 
jump instruction. 

We note that, generally speaking, an one-to-one correspondence exists between 
three foregoing techniques for describing computative algorithms. In other words 
these techniques are identical in substance. 

One of modern techniques available for describing computative algorithms is 
using the built-in predicates calculus, realised, for instance, in Prolog *. In particular, 
we may represent Prolog description of Euclidean algorithm by three statements 


GCD(0,V,V). (2) 
GCD(NS,VS,V):-NS1 is VS mod NS, GCD(NS1,NS,V). 
2-GCD(b, a, V). 


Where the second statement is the direct record of the recursive computative 
procedure, allowing to find GCD(a, 5); the first one determines the condition to 
finish this procedure; the third one is constructed to introduce the concrete values of 
numbers of Euclidean algorithm; NS, NS1, V and VS are internal variables of the 
procedure and V = GCD(a, b) after calculations. The main obstacle of this technique 
spreading is necessity of preliminary good knowledge of predicates calculus theory. 

This paper is devoted to an advance of analytical approach to describing some 
combinatorial and number-theoretic computative algorithms. Since at present any 
analytical description of the computative algorithms allows to automate the process 
of obtaining the program code, we suppose that the discussed theme appears to be 
interesting. 


2 Constructing analytical formulae by using logical functions 
2.1 Formulae to generate n-th prime number p, 


In our view, the most impressive application of logical functions in elementary 
number theory is the formula *’® to generate n-th prime number p,: 


(n+1)*+1 m 

Po=  san+1— Di (ED? ALE DY? / ED), @) 
m=0 k=2 

where po = 2, pi =3, ...; square brackets mean the integer part; sg is a logical 


function: sg(x)=1 if x>0 and sg(x)=0 if x <0. Let us find another analytical 
formula for p, without factonials. 
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As well-known” ’, any prime number has exactly two divisors: the unit and 
itself. Thus, any integer number a is a prime one if it has not divisors among integer 
numbers from 2 to [Va ] or, in the language of analytical formulae, if 


va 
) Thsee- sia )} (4) 
jr2 


then a is a prime only when x, = 1. It appears directly from (4) and (3) that the 
desirable formula for p, has form 


(n+1)7+1 m : 
P,= ¥ sgn-1- Dx); (5) 
33 


m=0 


where p2 = 2, p3 = 3, pa=5,.... 


2.2 The analytical description of the permutations generator 


As well-known”? 

a) the permutation of order m is called any arrangement of m various objects in a 
series; 

b) the verbal description of the simple algorithm available for constructing all the 
permutations from m objects, if all the permutations from m — 1 objects have been 
already constructed, has form 


Enumerate m-—1 various objects by the numbers 1, 2, .... m—1. For each 
permutation of a, a2, ..., a,_,, containing the numbers 1, 2, ...,m-—1, form m other 


permutations by putting the number m in all the possible places. As a result 
obtain the permutations: 


M, Qi, G2, ..5 An; (6) 
Qi, M, Q2, wy Anis 
GG GR Sve ran Sy oe Sao ; 
a, @2, .., M, as; 


a}, @2, .... @ 


It is evident that one can obtain all the permutations of order m by this algorithm 
and none of the permutations of (6) may be obtained more than once. If this verbal 
description we translate into one of special computer-oriented languages, for 
instance, into Pascal, then we obtain the program code, shown in table 1. This 
program works correctly at initial conditions m4=1; nl=m and the array nb3 
contains such numbers in the first m cells, which should be rearranged, and has 
following advantages over the verbal description 

a) the knowledge of all permutations from m — 1 objects is not required for 
generating the permutations of order m; 

b) permutations are realised with any set of numbers, contained in the first m 
cells of the array nb3. 
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Table 1. The representation of the permutations generator as Pascal program. 


Procedure Perm(Var m4,n1,n:Integer; If m2=n Then 
Var nb3,nb4,nb5:Ten); Begin 
Label A28,A29,A30; nb5{n]:=-1;Goto A29; 
Var nt,k,m2:Integer, End; 
Begin if Abs(m2)>0 Then Goto A30; 
lf m4=1 Then nb5[n]}:=1;ine(k); 
Begin A29: 
m4:=0;n:=n1; If n>2 Then 
For k:=2 to n do Begin 
Dec(n);Goto A28; 
End; 
Inc(m2);m4:=1; 
A30: 
m2:=m2+k; 
nt=nb3[m2]; 
nb3[m2]:=nb3[m2+1}; 
: nb3[m2+1]:=nt 
m2:=nb4[n}+nbS[n};nb4[n]:=m2; End; 


It should be noted that the main fault of both the verbal description and the 
program code is the fact that the knowledge of the previous permutation of order m 
is required for constructing the next permutation from m objects. To eliminate this 
fault one may use a set of analytical formulae ? 


n=pj-atl, p=j-1+f(l-g)+¢q(m-j-f), (7) 
S= tj (m-j+1) [t,,/(m -j+1)], 
bE e=g201 g=len4-1172, 
g=1 
z= sg(1+ p,_;— pj — 22) +22, 22= sg(1+ pj_.—pj-23) +23, «.- 
zj-1 = Sg(1+ pi — pj), 


where k is a number of permutation, generated of (7); 7; is a number, which j-th 
element of the initial sequence nb3 has in k-th permutation; the all another 
parameters in (7) are auxiliary. 


2.3 The formula for counting the value of GCD(a, b) 


We may present one of possible formulae available for counting the value of 
GCD(a, b) {see Sect. 1} as 


GCDca, b) = b{1 -sign(r)} + ksign(r), r=a-b[a/b}, (8) 


[vb] 
k= MAX fi(l-d)}, d=sign{a—ifa/i]} x sign{b — i[b/i]}, 
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where the function MAX(ai, a2, ...,a@,) gives the greatest from numbers a), a, ..., 
a,; sign(x) =|x\|/x if x 40 and sign(0) = 0. 


2.4 Formulae for the calculation of n-th numbers in Smarandache sequences of 
lst and 2nd kinds 


As we found earlier i the terms of six Smarandache sequences of Ist kind® may be 
computed by means of one general recurrent expression 

Gin) = (an 10¥) +a, +1), (9) 
where (nm) and (a,) are some functions; 6 is an operator. For instance, 


a) if p(n) = n+1, o=1 and y(a,) = [Ig(n + 1)] + 1 then the following sequence 
of the numbers, denoted as S;sequence, is generated by (9) 


1, 12, 123, 1234, 12345, 123456, ... (10) 
Let each number x; , determined as 
fig(k+0,5)] 
Xe =-l+ DL(k+1-10/), (11) 
j=0 


correspond to each number a; of sequences (10), where the notation “[lg(y)]” means 
integer part from decimal logarithm of y. Using (11) we may construct the following 
analytical formula for the calculation of m-th number in the S;-sequence: 


n 
a, = 10%" ¥(i/10%); (12) 
i=l 

b) if p(n) = n+1; o =yis the operator of mirror-symmetric extending the number 
41(n+1)/2] Of Si-sequence from the right with 1-truncating the reflected number from 
the left, if 7 is the odd number, and without truncating the reflected number, if n is 
the even number; y(a,) = [ lg({(m+1)/2] + 1)]+ 1, then the following sequence of 
the numbers, denoted as Sp-sequence, is generated by (9) 

1, 14, 121,-1221,. 12321;.123321, 1234321... (13) 


The analytical formula for the calculation of n-th number in the S2-sequence has the 
form 


[n/2] [(n+1)/2] 
= S10 “Hel + Siro? (14) 
i=1 i=l 


where d=1+ X[(n4+1)/2] + X[n/2] —Xi> and so on. 


We find recently that the terms of Smarandache sequences of 2nd kind® may be 
computed also by the universal analytical formula {compare with formulae (12) and 


(5)} 
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U, m 
an= > sg(n+2-b- > y;), (15) 
m=] i=b 
where %; are the characteristic numbers for Smarandache sequences of the 2nd 
kind; U, is an up-estimation for the value a,; b is a constant. For instance, if 


g! [Ve] 
Un=n’, b=2, x; =se{ >, []se(c-aic/g])}, g=flg+l, (16) 
k=1 q=l 


& 
c= 108 ¥ { [¥/108-? j-10 [4/108 ?*1} 3/107 , 
p=l 
the following sequence of Smarandache numbers of 2nd kind is generated by (15) 
1, 2,3, 5,7, 11, 13, 14, 16, 17, 19, 20, 23, 29, 30, 31, 32, 34, ... (17) 


We note that 

a) in formula (16): k is a number of the permutation, which is generated from 
digits of the number i; 7; is a number, which j-th digit of the number 7 has in k-th 
permutation {see (7)}; 

b) only such integer numbers belong to the sequence (17), which are prime 
numbers or can be derived to prime numbers by a permutation of digits in the initial 
natural numbers {the number 1 is related to prime numbers in this sequence}. 


2.5 Formulae for analytical description of Magic squares constructing methods 


As we discovered earlier® '', logical functions may be used effectively to create 
analytical formulae available for describing computative algorithms on constructing 
Magic squares of any order n. For instance, let us consider a well-known “Method of 
two squares”, whose verbal description has the form*”*: 


1. Make a drawing of two square tables of any order n = 4k +2 (k=1, 2, ...). 
Divide every table in four equal squares which we shall call A, B, C and D squares 
respectively {see figure 2(1)}; 

2. Place a Magic square of order m= 2k +1 inthe A, B, C and D squares of the 
first table. It is obvious {see figure 2(2)}, the first table will have the same sum of 
numbers in its rows, columns and main diagonals; 

3. Fill the cells of the second table: all cells of A square are to have zeros; cells 
of D square — numbers u = m’; cells of B square — numbers 2u and cells of C 
square — numbers 3u. The obtained table {see figure 2(3)} will have the same sum 
of numbers only in its columns; 

4. Perform such rearrangement of the numbers in the table 2(3) that the new table 
will have the same sum of numbers in its rows, columns and main diagonals. It can 
be achieved, for instance, by operations 

a) underline in the square A of the second table 

— k zeros, located in the extreme left positions of all rows, excepting the middle 
row {see figure 2(3)}; 
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— the zero, located in the central cell of the middle row, and another k— 1 zeros, 
located left of the central cell. 

Exchange all marked zeros against the respective numbers of the square C {see 
figure 2(4)} and otherwise; 

b) mark k-1 numbers 2u in the extreme right positions of every row of 
square B {see figure 2(3)}, and then exchange ones against corresponding numbers 
of the square D (see figure 2(4)) and otherwise. 

5. Add (cell-wise) two auxiliary tables {the Magic square of order 10, obtained 
as aresult of adding auxiliary squares 2(2) and 2(4) is shown in figure 2(5)}. 


75 | 75 | 75 | 75 | 75] 
75 | 75 | 75 | 75 | 75 
75 | 75 | 75 | 75 | 75] 


25 | 25 | 25 | 25 | 25 
25 | 25 | 25 | 25 | 25 | 
75 | 75 | 75 | 75 | 75] 25 | 25 | 25 | 25 | 25 | 
75 [75 | 75 | 75 | 75} 25 | 25 | 25 | 25 | 25 


(3) 


71 | 52 | 68 | 60 | 34| 
3 | 55 | 74 | 51 | 42 | 
}7a{91 | 12| @ | 201] 54 | 66 | 62 | sa | 45 
j22 98 | 0 | 19 11] 57 | 73 | so | 6 | 6 | 
[+514 [01 | 97| ral] 40| 39 | 91 | 47 | 59 | 
j21| 2 | sa | as | ofa | 27 | 49 | 35 | 59 
jee | 5 | 24 | 76 | sal] 30 | 20 | 49 | 26 | 67 
[4 | 16] 67 | 63] 9529 | 41 | 97 | 29 | 70 | 
[0 | 0} 75] 75] 75H} 2s | 25] 25] 25] 50] | 7 | 29] 75 | 94 | oo] 2 | 40 | 25 | «4 | | 
(4) (5) 


Figure 2. Method of two squares. 


[0 | 0 | 75] 75 75| 
jo | o | 75| 75 | 75 
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Table 2. The representation of computative algorithm “Method of two squares” as Pascal program. 


Procedure BuildMS; {Modify the auxiliary square } 
{ Build Magic square of order n=4k+2 } { Modify A and C squares } 
Var m,i,j,n2:Byte; For i:=1 to n2 do 
zZ,.x,y,nn:integer; tf i<>n2 Shr 1+1 Then 
Begin Begin 
2:=0;n2:=n; For j:=1 to m do 
n:=n2 Shr 1; Begin 
m:=n Shr 1; Inc(mk{Ind(i,j)],an*3); 
nn:=Sar(n); Dec(mk[Ind(i+n2,j)],nn*3); 
{ Build an auxiliary square of order n=4k+2} End; 
{ Place a Magic square in A square } End; 
For z:=1 to nn do i:=n2 Shr 1+1; 
Begin For j:=i Downto i-m+1 do 
x:=(fX(z}-1+n)Mod n+1; Begin 
y.=(gY(z}-1+n)Mod n+1; Inc(mk{Ind(i,j)],nn*3); 
mk[(x-1)*n2+y]:=z; Dec(mk[Ind(i+n2,j)],nn*3); 
End; End; 
n:=n2;n2:=n Shr 4; { Modify B and D squares } 
For i:=1 to n2 do For j:=1 to n2 do For j:=n Downto n-m+2 do 
Begin For i:=1 to n2 do 
{ Place the Magic square in B, C and D squares } Begin 
mk({Ind(i,j+n2)}:=mk{Ind(i,j)}+nn Shi 1; Dec(mk{Ind(i,j)],nn); 
mk{Ind(i+n2,))]:=mk{Ind(i,})]}+An Shi 1+nn; inc{mk{Ind(i+n2,j)],nNn); 
mk{Ind(i+n2,j+n2)}:=mk{ind(i,j)}+nn; End; 
End; End; 


If the foregoing verbal description we translate into Pascal, we obtain the 
program code, shown in table 2. In this program code 

a) the Magic square of order m = 2k + 1, located in the A, B, C and D squares of 
the auxiliary square, may be built for instance, by the functions »'> 


fiz, m) = (z-1) + [@-1)/m] - [m/2], (18) 
g(z, m) = (¢-1)- [(@-1)/m] +[m/2], 


where square brackets mean the integer part; a sign “=” is the modulo m equality; z 
is any natural number from 1 to m’; functions f and g afford to compute the position 
of any natural number z in cells of the Magic square: x = f(z, n) and y = g(z, n). In 
particular, functions (18) may be presented as following two distinct Pascal- 
procedures 


Function fX(z:integer):Integer, Function gY(z:Integer): integer; 
Begin Begin 

fX:=1+(z-1)+(z-1}div n - n shr 1; gY:=1+(Z-1 {z-1)Div n +n shr 1; 
End; End; 


b) two procedures “Ind” and “Sign” are auxiliary and have the form 


Function Ind(x,y:Integer):integer, Function Sign(n:Word):Shortint; 
Begin Begin 

Ind:=(x-1)*n+y; tf Odd(n) Then Sign:= -1 Else Sign:=1; 
End; End; 
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The analytical description of “Method of two squares” has the form? 


x =item—-1; y=(l—-c3-—cg—cs) (jt+e.m) + (e3+c6+es) (1+ (19) 


+ {7+ (c2 + Im — 1) mod n}) - 1; 


u=m’; z=1+ {(W-1)modu}; i=f(z, m) + 1; 


J= gZm)+1, cr =[{([V-1)/u) + 1) mod 4}/2 }; 


C2 


[(V-1)/u]mod2, c3=[ (sign(c} m+ i-3k-4) + 2)/2 ]; 


ca = asg(j—k-1); ceo=c4[ (sign(k-—i-cym) + 2)/2 }; 
5 = (1 —c4)(1 —c1)[ {sign — 1}+ 1}/2 Jx[ (sign(k-7 + 1) + 2)/2], 


where n = 4k + 2 is an order of the desirable Magic square, contained natural 
numbers N from | to n’; m= 2k + 1; functions f(z, m) and g(z, m) are determined 
by the expression (18); asg(x)=1 if x#0 and asg(0)=0 {asg(x) = |sign(x) | 


= $1 


X i= SHxX |G. 
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